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Abstract 

We investigate the isotropic-anisotropic phase transition of the two- 
dimensional XY model with six-fold anisotropy, using Monte Carlo renormal- 
ization group method. The result indicates difficulty of observing asymptotic 
critical behavior in Monte Carlo simulations, owing to the marginal flow at 
the fixed point. 
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Typeset using REVTgX 



The XY model with Zq symmetry breaking field represents not only planar spin magnets 
with Z 6 symmetric crystal field but also models with Z 2 X Z 3 symmetry (which is 
isomorphic to Zq) such as three-state antiferromagnetic Potts model on the square and 
the cubic lattice and Ising antiferromagnet on the triangular and the stacked-triangular 
lattice g. 

The critical behavior of the model in two dimension is well understood by scaling 
argument [|TJ: The model undergoes two distinct phase transitions, both of which be- 
ing Kosterlitz-Thouless transition. In Monte Carlo simulations of finite system, however, 
marginal renormalization flow near the fixed point makes it difficult to observe asymptotic 
critical behavior. In the present work, we used improved Monte Carlo renormalization group 
method |§ to observe the renormalization flow of the six-clock model on the square lattice 
of size up to 64 x 64, and found that the size 64 x 64 is still insufficient to observe asymp- 
totic critical behavior. Similar situation occurs in several other models, such as general spin 
models in D = 4 || and 0(3) spin model with cubic anisotropy in D = 3 [IT]]. 

In the Ref. |J by Jose et.al., they investigated the following model: 

H = KY. cos(9i - 6 3 ) + cos( P 6i) (1) 

<ij> i 

where the first summation runs over all nearest pairs on the square lattice, and p is some 
integer. From Gaussian spin-wave theory and scaling analysis, it can be shown that the 
perturbation term J2iCos(p8i) is relevant when rj < 4/p 2 , where t] is the critical exponent of 
long distance spin-spin correlation defined as below: 

< cos(8 x - 6y) >~ \x - yl^ 71 for |x-y|>l. (2) 

Thus flow of the renormalized parameters K and A for p = 6 case is expected to become 
the one shown in Fig. [I]: The parameter space is divided into three regions, namely high- 
temperature (H), Kosterlitz-Thouless (KT), and low-temperature (L) phases. Any model 
whose initial parameter crosses the L-KT boundary is attracted to a fixed point which is 
on the KT fixed line and characterized by an exponent r\ = 1/9 (we denote the fixed point 
as F\/g) and exhibits KT-phase transition. Near the fixed point i^i/g, renormalization flow 
along the L-KT boundary line slows down because Fim is marginal for both direction K and 
A. Thus finite size correction (distance to the final fixed point) is expected to behave like 
1 / log L (L is linear size of the system) and one can not observe asymptotic critical behavior 
unless extremely large system is used. For example, the critical exponent rj at the lower 
critical temperature of Ising antiferromagnet on the triangular lattice has been estimated by 
several authors (to confirm theoretical prediction v = 1/9 = 1.111 • • •) as r\ — 0.15(2) Q and 
r\ = 0.125(25) 0. The accuracy of these values are relatively low compared to other models 
such as 0(n) spin models in three dimension: this implies that there are large finite-size 
correction. 

In the present work we numerically investigated the renormalization flow of the model 
(|I]) using the improved Monte Carlo Renormalization Group (MCRG) method 0, which 
is very simple and efficient way to extract essential information of critical phenomena from 
simulation data. We observed the following quantities: 

< S(kx) ■ S(-kQ > <i&cos(60M)> ^ 

KL=1 ^ ' < >3 (3) 
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where 



S(k) = L~ 2 5>xp(ik-x)(cos# x , sin# x ), ki = (2tt/L, 0), (4) 

X 

and 

M = S(0) =R M (cos6 M , sm9 M ). (5) 

Kl and Al reflect the behavior of the renormalized temperature and anisotropy (by a factor 
L), respectively. We also observed Binder's parameter and < cos 69m >> an d found that K L 
and Al reflect the RG flow better than these quantities. 

The Hamiltonian (|) on a L x L square lattice is simulated for L = 16, 32, and 
64. In Monte Carlo simulation, Metropolis update scheme and Wolff's cluster algorithm 
]nj are combined. In the Metropolis update, we choose new spin value 9 with a prob- 
ability proportional to exp[— A cos(60)] and calculate acceptance probability using only 
KJ2<ij> cos(9i — 9j). This scheme satisfies the detailed-balance condition of the Hamiltonian 
(H) and improves the acceptance ratio for large A. In the cluster update, spin-reflection axis 
is restricted to the one which preserves the anisotropy. The Hamiltonian (P is simulated 
for A = oo, 0.5, 0.2, 0.1, and 0.05 cases. 

Figure || shows the flow of K L and A L . each line is drawn from (K L , A L ) to (K 2 l, A 2 l)- 
One can see that the obtained RG flow agrees well with the theoretical one. The position 
of the critical point Fi/g is estimated by plotting log < M| > +(logL)/9 versus logL for 
various K, being A fixed to zero: since < M| >~ L~ v , the plot becomes horizontal at -F1/9. 
Fig. |^ indicates that F1/9 is located near K = 1.7. The value of Kl at this temperature 
is marked as Fi/g in Fig. [2]. One can see that the RG flow is attracted to the KT line in 
the left (high-temperature) side of F1/9, while the flow deviates from the KT line in the 
right (low-temperature) side of F1/9. Plots for A = 00 case seems to approach F1/9 as L 
increases. However, the approach is indeed slow and extremely large L is required to observe 
convergence to F 1 / 9 . 

The distance to F\/g results in systematic error in the finite-size scaling analysis [|],[|, 
such as position of the critical point and value of the critical exponent, and one can not 
extrapolate the L = 00 limit owing to the slow vanishing correction term. For example, if 
one define finite-size critical point as a temperature where renormalized anisotropy becomes 
size-independent, it will deviate to low-temperature side as Fig. [I] suggests. Similarly, the 
1] = 1/9 criterion will lead to systematic error, because effective exponent for rj at A ^ 
region generally differ from that of asymptotic value at A = 0. 

In conclusion, we demonstrated that it is nearly impossible, by means of Monte Carlo 
simulations of finite system to observe asymptotic critical behavior of isotropic-anisotropic 
phase transition in 2D six-clock model and those of similar symmetry. Although quantitative 
information such as position of critical point or values of critical exponents are hard to obtain, 
qualitative information such as presence or absence of the transition can be easily obtained 
by observing renormalized anisotropy, and it is enough for Monte Carlo simulation since the 
nature of the transition, such as critical exponents, is already well known. 
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FIGURES 
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FIG. 1. Theoretical prediction of the renormalization flow of temperature and anisotropy. L, 
KT, and H denotes low-temperature, Kosterlitz-Thouless, and high-temperature phase, respec- 
tively. Dotted lines indicate direction of marginal renormalization flow. 
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FIG. 2. Renormalization flow obtained from Monte Carlo simulations. Connected two line 
segments are drawn from L = 16 data to L = 32 data, then to L = 64 data. Other single lines 
are drawn from L = 16 data to L = 32 data. Plots for each different values of A are separated by 
dotted lines. The position of F 1 / 9 was estimated from finite-size scaling of < M\ >. 
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